SPECTRAL DEVIATIONS FOR THE DAMPED WAVE EQUATION 
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Abstract. We prove a Weyl-type fractal upper bound for the spectrum of the damped 
wave equation, on a negatively curved compact manifold. It is known that most of the 
eigenvalues have an imaginary part close to the average of the damping function. We 
count the number of eigenvalues in a given horizontal strip deviating from this typical 
behaviour; the exponent that appears naturally is the 'entropy' that gives the deviation 
rate from the Birkhoff ergodic theorem for the geodesic flow. A Weyl-type lower bound 
is still far from reach; but in the particular case of arithmetic surfaces, and for a strong 
enough damping, we can use the trace formula to prove a result going in this direction. 



1. Results and questions about the spectrum of the damped wave 

equation 

Let (M, g) be a smooth compact Riemannian manifold without boundary. Let a be a 
C°° real valued function on M. We study the "damped^ wave equation", 

(1.1) (dl - A + 2a(x)d t ) v = 

for t G R and x G M. We shall be interested in the stationary solutions, that is, solutions 
of the form v(t,x) = e ltT u(x) for some r G C. This means that u must satisfy 

(1.2) (— A — r 2 + 2iar)u = 0. 
Equivalently, r is an eigenvalue of the operator 

I 

-A 2ia 

acting on H l (M) x L 2 (M). For a = this operator is the wave operator, an anti-selfadjoint 
operator; but for a ^ the operator is not normal anymore. It is known that its spectrum 
is discrete and consists of a sequence (r n ) with 3m(r„) bounded and |3?e(r n )| — ► +oo (see 
Section [2|). One sees easily that Qm(r n ) G [2 min(inf a, 0), 2 max(sup a, 0)] if 9fte(r n ) = 0, 
and ^sm(T n ) G [inf a, sup a] if 3?e(r n ) ^ [28]. One can also note, obviously, that — f n is an 
eigenvalue if r n is : the spectrum is symmetric with respect to the imaginary axis. 



The term "damped" applies to the case when a > 0, that is to say, when the energy is decreasing. In 
this paper the sign of a will be of no importance. 
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1.1. Background. Motivated by questions from control theory, Lebeau [17] was interested 
in the so-called "stabilization problem" : define 

p = sup {(3,3C, E{u t ) < Ce- 0t E(u o ) for every solution u of (Oi l} 

where the energy functional is E{u) = f M \ Vu\ 2 . The stabilization problem consists in 
finding some damping function a (necessarily nonnegative) such that p > 0. Lebeau 
identified 

p = 2 min {£>(()), C(oo)} 
where D(0) is a sort of "spectral gap" : 

D(0) =inf{3m(r n ),T n ^0} 
and C(oo) is defined in terms of the Birkhoff averages of a along the geodesic flow 

C : T*M — ► I'M. 

-t 



i r 

C(oo) = lim inf - / a(G s p)ds. 

t — >+oo per* m t J 



Lebeau also showed, on an example, that it is possible to have a spectral gap (-D(O) > 0) 
but no exponential damping (p = 0) : this surprising phenomenon is typical of non normal 
operators. 

Markus and Matsaev [19] proved an analogue of Weyl's law, also found independently 
later on by Sjostrand [28] : 

(1.3) Un,0<Xe(T n )<\}=(£-) ([ dxdi + 0{\- 

where d is the dimension of M, p is the principal symbol of —A, namely the function 
p(x, £) = g x (£,Q defined on the cotangent bundle T*M, and dxd^ is the Liouville measure 
on T*M (coming from its canonical symplectic structure). 

It is a natural question to ask about the distribution of the imaginary parts 3m(r n ) 
in the interval [inf a, sup a]. For non normal operators, obtaining fine information on the 
distribution of the spectrum is much harder than for normal operators - one of the reason 
being the absence of continuous (or even smooth) functional calculus. Another related 
difficulty is that there is no general relation between the norm of the resolvent and the 
distance to the spectrum. Some techniques have being developed to obtain upper bounds 
on the density of eigenvalues, but lower bounds are notoriously more difficult. 

Assuming that the geodesic flow is ergodic with respect to the Liouville measure on 
an energy layer, Sjostrand proved that most of the r n have asymptotically an imagi- 
nary part 3m(r n ) ~ a, where a denotes the average of a on the energy layer S*M = 
{(x, £) G T*M, g x (£, £) = 1} with respect to the Liouville measure : 

Theorem 1.1. [28] Assume that the geodesic flow is ergodic with respect to the Liouville 
measure on S*M. For every C > 0, for every e > 0, we have 

Un,X< ^e(r„) < A + C, 9m(r„) £ [a - e, a + e]} = o^ 1 ) 

as X goes to infinity. 
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Remark 1.2. If C is not too small, one sees from fl 1 . 3H that there exists C > such that 

H {n, A < 3?e(r n ) < A + C} > CA^ 1 
for large A. Thus, the theorem does say that a majority of the r n have 3m(r n ) G [a— e, a+e]. 



Remark 1.3. More generally, without the ergodicity assumption, [28j proves the following 
results. Introduce the functions on T*M, 

"T/2 



1 P 7 

(a) T = 7F / a ° GS(is ' 

J 7 -T/2 

(a)^ = lim (a) r , 



and the real numbers 



a + = lim sup (aW, 
a_ = lim inf (aW. 

T — ►+ooS*Af 

The function (a)^ is well defined Liouville-almost-everywhere, by the Birkhoff theorem. 

Lebeau [TTj proves that for any e > 0, there are at most finitely many n with r n $ 
[a_ — e, a + + e], and Sjostrand [28J proves that 

jj {n, A < 9fte(r n ) < A + C, 3m(r n ) G" [ess inf (a)oo — e, ess sup(a)oo + e]} = o(A -1 ). 

1.2. Semiclassical formulation. As in [28] we reformulate the problem using semiclas- 
sical notations. In doing so, we also generalize a little the problem. We introduce a 
semiclassical parameter < h <C 1 and will be interested in the eigenvalues r such that 

tvr — > 1. If we put r = ~ with A close to 1, then equation (11.21) can be rewritten as 

r — >o v ' 

/ fc2^ \2 \ 

(1.4) f — - — + ih\a)u = 0, 

or 

(1.5) (V-z)u = Q 
if we put 2 = ^ , and 

ft 2 A 

(1.6) V = V{z) = P + ihQ{z), P = — , Q{z) = a^. 

The parameter z is close to E — |. 

More generally, we will consider a spectral problem of the form 1 . 51) where 

7> = V{z) = P + ihQ(z), P = 

z e n = e 4 ]-*°' s °[]£ min ,£ max [, with < E min < \ < £ max < +00, < s < f. We will 
assume that Q(z) G ^DO 1 is a pseudodifferential operator that depends holomorphically 
on z G Q, and that Q is formally self-adjoint for z real (the definition of our operator 
classes is given in Section [9j). 
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We denote 

£ = Z h = { z G Q, 3u, (P(z) - z)u = 0} . 

The operator P has principal symbol p (x, £) = 9 *(|'^ , and Q(z) has principal symbol 
taking real values for 2 real. In these notations, the previous results read as 
follows : for any E min < E 1 < E 2 < E mSuX , 



(1.7) ft i z e S, Ex < Ke(z) < E 2 } ' 



(27rft) £ 



+ 0(H) 



(uniformly for all c > and for all Ei,E 2 such that \E 2 — Ei\ > 2ch, E\ and E 2 staying 
away from E min , E m3iX ). One can show that (for z G £) has to stay bounded if £a, i?2 

stay in a bounded interval. Fix some c > 0, and take E\ = E — ch and E 2 = E + ch. Let 
us denote 

Qe = „ lim i nf (Qe)t, 

9£ = lim sup (g E ) r , 
r -" +0O Po - 1 {s} 

then we have [28] 

(1.8) g - + o( l)<__^< g + + (l) 

for all z G £ such that E — ch < 9fte(z) < E + ch. Finally, denote q~E the average of on 
the energy layer p^ 1 {E}. Assuming that the geodesic flow is ergodic on p~ 1 {E}, then for 
any e > 0, any c > 0, 

(1.9) ft jz G £, £ - eft < 5fte(z) <E + ch, £ - e, + e] j = o^). 

The method of [28j allows to treat the case of a more general P (and thus a more general 
Hamiltonian flow than the geodesic flow), and also to deal with the case when the flow is 
not ergodic. However, in this paper we will stick to the case of an ergodic geodesic flow; 
we will add even stronger assumptions in the next paragraph. 

1.3. Questions, and a few results. We try to give (partial) answers to the three natural 
questions : 

(Ql) (Fractal Weyl law) Let / be an interval that does not contain qg. Can we 
describe in a finer way the asymptotic behaviour for 

z,E — ch< $te{z) <E + ch, ^iM. e / J ? 

For instance, can we find 

logftjz G £, E - ch < fte(z) <E + ch,^l E l\ 

lim — ? 

h cQ log h 
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(Q2) (Inverse problem) Suppose we know everything about P and about the dynamics 
of the geodesic flow, but that Q is unknown. To what extent does the knowledge of the 
eigenfrequencies {z G E, E — ch < §te(z) < E + ch} determine qE ? 

Replacing V by B~ X VB^ where B is an elliptic pseudodifferential operator with positive 
principal symbol 6, amounts to replacing q by q — {p Q ,\ogb}, where {.,.} stands for the 
Poisson bracket on T*M. Two smooth functions / and g on T*M are said to be cohomolo- 
gous (with respect to the geodesic flow) if there exists a third smooth function h such that 
/ = g + {p , h}. This defines an equivalence relation, we write / ~ Po g. 

It is thus more natural to ask : 

(Q2') Does the knowledge of the eigenfrequencies {z G E, E — ch < %te{z) < E + ch} 
determine the cohomology class of qE ? 

If the length spectrum of M is simple, then one can most probably use a trace formula 
and recover from the knowledge of E all the integrals of qE along closed geodesies. And this 
is enough to determine the cohomology class of qE if M has negative sectional curvature : 
the Livshitz theorem [HI [7] says that if two functions have the same integrals along all 
closed geodesies, then they are cohomologous. Thus, the answer to (Q2) is probably "yes" if 
M has negative sectional curvature and the length spectrum is simple (this last assumption 
is satisfied generically, but unfortunately not for surfaces of constant negative curvature) . 

In fact, it also makes sense to ask whether some knowledge of the imaginary parts alone 
{3m(z), z G E, E — ch < ?fce(z) < E + ch} allows to recover some information about qE- 
For instance, one can ask the apparently simple question : 

(Q3) (Very weak inverse problem) Let C denote a constant function. As follows from 
(11.81) . we have the implication 



Does the converse hold ? 

Main assumption on the manifold M : From now on, we assume that M has 
constant sectional curvature —1. This implies the ergodicity of the geodesic flow on any 
energy layer (with respect to the Liouville measure), and in fact a very chaotic behaviour 
of trajectories (see Section [3|). We will indicate how to generalize our results in the case 
of surfaces of variable negative curvature; however, it is not clear what to do in higher 
dimension and variable negative curvature. 

In the following theorem, hxs stands for the Kolmogorov-Sinai entropy, or metric en- 
tropy. It is an affine functional, taking nonnegative values, defined on Ai, the set of 
G-invariant probability measures on T*M : see for instance [12] for the definition of 
hxs- We will also denote M.e C M. the set of invariant probability measures carried by 
p~ 1 {E}. Since p Q is homogeneous it is enough to consider, for instance, the case E — =, 
and p~ l {E} = S*M. For /i G Ml, we have hxsifi) 5; d — 1, with equality if and only if 




C. 
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\i is the Liouville measure. We now fix E = | and we denote q — qi , q — q~i , q + — qt 

2 2 2 

<r = Qi- 

2 

We fix some c > and denote 

Ei = lz E S, ~ - ch < $te(z) < ~ + ch\ . 

Theorem 1.4. Assume M has constant sectional curvature — 1. Define 

H(a) = sup fi E -M-ii j qdjj, = a^. 

If ex > q, then for any c > 



limsup FT~Tfc1 — H( a )- 

h — 

If ot < q, then for any c > 





■ {zeEi 


Sm(z) 

' h 


> «} 




|log 


h\ 






: |z G Si 


^sm(z) 
' ft 


< a} 




|log 


h\ 





limsup — — — — < H(a). 

h — >o 

Remark 1.5. (see [15] . §4, or [4], §3 for the argument) The function H is concave and 
is identically — oo outside [q~,q + ]. It is continuous and strictly concave in [q~,q + ], and 
real analytic in ]q~, q + [ (note that q~ = q + if and only if q is cohomologous to a constant 
on S*M). The function H reaches its maximum d — 1 at the point g, and has finite 
(nonnegative) limits at the endpoints q~,q + . In particular H is positive in the open 
interval }q~ , q + [. 

Remark 1.6. The key fact in the proof of Theorem 11.41 is the large deviation estimate 
from [T3] , 

log Li {p E S*M, (o)t{p) £ 1} (TTf n j 

^ lim — = sup {H(a), a E I \ — a — 1 

for any interval / C R - where Li is the Liouville measure on p~ 1 {|} = S*M. This result 
gives the volume of the set of trajectories deviating from the Birkhoff ergodic theorem. See 
Section [31 

On a surface of variable negative curvature, we can generalize the result to : 

Theorem 1.7. Assume M is a surface of variable negative curvature. Denote tp the 
infinitesimal unstable Jacobian (see Section^). Define 

ui \ ( h K s{v) ^ . A f , 
H (a ) = sup < -= -,/jGMi, / qd[i = a 

U (fdfi 2 J 

Then the same conclusion as in Theorem \l.J\ holds. 
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Remark 1.8. For a manifold of variable negative curvature and dimension d we would 
expect the same to hold with H(a) = (d — 1) sup j j^f ^ , /i€ MiJ g <i/i = a j . However, 
our proof does not work in this case. 

One may wonder if the limsup in Theorem 11.41 is also a liminf. This question is rem- 
iniscent of the fractal Weyl law conjecture asked by Zworski and Sjostrand, but in our 
situation we can say with certainty that it is not true. Worse, one cannot expect the lower 
bound to hold for a "generic" q. In a paper in preparation, Emmanuel SchenckH is proving 
that there exists 5 > such that < q + — § for all z G Si, provided a certain criterion 

on q is satisfied. The criterion reads Pr(g — < g + , where the pressure functional Pr 
is defined on the space of continuous functions on S*M and is the Legendre transform of 
~hxs ( see Section [3]). The functional Pr is lipschitz with respect to the C° norm, and the 
condition Pr (q — < q + defines a non-empty open set in the C° topology. For such a 
q we cannot have 



logtf{ze£i,^P >q + -d) 



liminf — > H(q+ - 5), 

h^o \hgh\ 

since H is positive in ]g_,g + [ but the liminf on the left-hand side is — oo. 

In Sections [7] and El we investigate Question 3 in some special cases. We work on compact 
hyperbolic surfaces (d — 2), and we study operators of the form 

A u f = Af-2(u,df)+\\u\\ 2 x f, 

called "twisted laplacians" - here to is a harmonic real-valued 1-form on M. Studying the 
large eigenvalues of A w amounts to studying a fixed spectral window for the semiclassical 
twisted laplacian 

_£2^ = - ft 2 — + £2 ( ±) _ h 2 JMg fa > q 

2 2 \ ' / 2 

This question falls exactly into the setting of §1.21 with q(x, £) = (u) x ,£). Since q(x, — ^) = 
— q(x, £), we have q = 0. We will denote Pr(cj) = Pr(g) the pressure of the function q, 
defined in Section I3.2L It will also be interesting to note that q + = —q~ coincides with 
the stable norm ||o;|| s defined in Section 17.21 that ||u;|| s + 1 > Px(uj) > \\u)\\ s and that 
Pr(tu) — \t\\\u\\ a — >■ in the case of surfaces. 

t >oo 

Theorem 1.9. Assume M is a compact arithmetic surface coming from a quaternion 
algebra. Take w ^ 0. Fix (3 G (0,1], and < e < (3. Then, for every h small enough, there 
exists z G Sp(-h 2 ^) with \Ue(z) - || < h 1 ' 13 , such that 

^sm(z) s 1 1 + e 

> Pr (uj 



h ~ v ' 2 2(5 



2 Private communication. 
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Equivalently, given (3 £ (0, 1], and < e < (3, for all T large enough, there exists r n such 
that |»e(r n ) — T\<T 13 and 

\Qm(r n )\ > Pt(uj) - - - 

where r n is the "spectral parameter" defined by X n = — (| + r%). 

Of course, we deduce immediately the following corollary : 

Corollary 1. H{n, |3?e(r n )| < T, |Qm(r n )| > Pr(cj) — | — ^} > T 1_/3 , asymptotically as 
T — ► +oo. 

Since Pr(u) — > 1 as uj — ► (but Pr(cj) — ► +oo as ||a;|| s — ► +oo) this result is only 
interesting if ||a>|| s is large enough (depending on (3). Note also that if ||a>|| s is large enough 
we have Pr(a>) < ||a;|| s + |, so that the work of Emmanuel Schenck mentioned above will 
show that there is a strip (Sm(z) > \\uj\\ s — 5} (5 > 0) that contains no r n . 

The arithmetic case is special, in that the lengths of closed geodesies are well separated : 
we can write a trace formula, find a lower bound on the geometric part (despite of its 
oscillatory nature) and deduce a lower bound on the imaginary parts of eigenvalues. The 
ideas are borrowed from |8j. 

Remark 1.10. Another way of defining the twisted laplacian is to write M = T\M, where 
EI is the universal cover of M and T is a discrete group of isometries of H; to fix an origin 
o £ H, and to write 

AJ(x)= e />oA(e^>/W); 

this operator preserves T-periodic functions, hence descends to M. The case where to takes 
purely imaginary values is self-adjoint, and analogous to the study of "Bloch waves". The 
case where ui takes real values is no longer self-adjoint. 

Remark 1.11. Our motivation for working with twisted laplacians on hyperbolic manifolds 
was that it is a case where the trace formula is exact. There was, a priori, hope to prove 
finer results than in cases where the trace formula is not exact (in the latter case the space 
of test functions to which the formula can be applied is more limited). A posteriori, we 
never use any wild test function that wouldn't be allowed in the general case. So, one can 
think that the same result holds for our general operator (JTTHJ) - provided one proves a 
semiclassical trace formula first. 

If we don't assume arithmeticity, we can only treat the following operator : 

where 9(h) > \ \ogh\ and h9(h) — > 0. In other words the non-selfadjoint perturbation is 
stronger than in the previous cases. 
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Theorem 1.12. Assume M is a compact hyperbolic surface. Take uj 7^ 0. Take any 
function f(h) ^> 9(h) 3 / 2 log log h 1 ^ 2 . Then there is a sequence h n — > such that 

the stable norm of u. 

See Section 17.21 for the definition of the stable norm. Note that with our previous 
notations, ||u;|| s = q + = —q~, and Pr(iu;) — |t|||a;|| s — > on a hyperbolic surface. 

t >oo 

Acknowledgements : This work was partially supported by the grant ANR-05-JCJC- 
0107-01. I am grateful to UC Berkeley and the Miller Institute for their hospitaly in the 
spring of 2009. I thank Dima Jakobson for suggesting that some of the ideas contained 
in Hejhal's book [8] could be used to find lower bounds on the density of eigenvalues in 
arithmetic situations. In fact, at the same time as this paper was written, Jakobson and 
Naud managed to apply these ideas to study the resonances of certain arithmetic convex 
cocompact surfaces [TO] . 

2. Note on the definition of the spectrum and its multiplicity. 

By the term "spectrum", we mean the set £ of z G Q such that Viz) — z is not bijective 
H 2 (M) — > H°(M). If it is bijective, then the inverse must be continuous, by the closed 
graph theorem. 

If z is restricted to a compact subset of fl, it is easy to see that G(z) = I + X~ 1 (V(z) — z) 
is invertible for A > large enough. The inverse G(z)~ x is then a compact operator 
on H°(M). Besides, one sees that V(z) — z is not bijective H 2 (M) — ► H°(M) if and 
only if 1 is in the spectrum of G(z)~ l , if and only if there exists u G H 2 (M) such that 
(V(z) — z)u = 0. This shows, in particular, that the "spectrum" is discrete and corresponds 
to the existence of "eigenfunctions". 

To define the multiplicity of z G E, we proceed the same way as in [28]. By the 
density of finite rank operators in the space of compact operators [24], one shows that in a 
neighbourhood of z there exists a finite rank operator K(z), depending holomorphically 
on z, such that V(z ) — z + K{zq) is invertible. The multiplicity of z is then defined as 
the order of Zq as a zero of the holomorphic function 

z ^ det[(V(z) -z + K(z))-\V(z) - z)\ = det[J - (V(z) -z + K(z))- 1 K(z)}. 

It is shown in [28] that this definition does not depend on the choice of K(z), and coincides 
with other usual definitions of the multiplicity. Besides, the argument can be extended to 
the case where K(z) is trace class. 

3. A FEW FACTS ON THE GEODESIC FLOW ON A NEGATIVELY CURVED MANIFOLD 



3.1. Anosov property. If M has negative sectional curvature, then the geodesic flow 
on S*M has the Anosov property [3]. This means there are C, A > such that for each 
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p G S*M, the tangent space T p (S*M) splits into 

T p (S*M) = E u {p) © E s {p) © RX(p) 

where 

- the vector field X generates the geodesic flow G tm , 

- E s is the stable subspace : for all v G E s (p), and for t > 0, H-DGt- v\\ < Ce~ xt \\v\\; 

- E u is the unstable subspace : for all v G E u (p), and for t < 0, ||IXt* i> || < Ce A *||v||. 

If M has constant negative curvature —1, any A < 1 will do. We take A = 1 — e, with e 
arbitrarily small. One also has an upper bound ||-DCr*. u|| < Ce^ 1+e ^ \\v || for any e > and 
any tel. 

3.2. Pressure, entropy, and large deviation. The pressure is defined on C°(S*M), as 

the Legendre transform of the entropy : 



Pr(/) = sup j/ii^O) + J fdn,neM h \ 



If / is Holder, then the supremum is attained for a unique fj,, called the equilibrium measure 
of /. The functional Pr is analytic on any Banach space of sufficiently regular functions 
- for instance, a space of Holder functions [HI [26]. Besides, the restriction of Pr to any 
line {/ + tg,t G M} is strictly convex, unless g is cohomologous to a constant [23]. If g 
is sufficiently smooth, we recall that this means that g = {p Ql h} + c for some smooth 
function h and a constant c. If g if Holder, it is better to use the integral version of the 
notion. If j(t) is a periodic trajectory of the geodesic flow on S*M (equivalently, a closed 
geodesic), we denote Z 7 its period (equivalently, the length of the closed geodesic). We 
denote the measure J g dj = g(^{t))dt on S*M, and d/x 7 the probability measure 
J g dfiry = ly 1 J gd-f. One says that g is cohomologous to the constant function c if 
f g dj = for all periodic trajectories of the geodesic flow (the Livschitz theorem says 
that both notions are equivalent for smooth functions). 

Let us now fix a smooth function q on S*M, not cohomologous to a constant. For «6l, 
define 



H(a) = sup |/ir:s(aO) A* e J l^n 



a 

P{(3) = Pr(/3g) = sup ^h KS (p,) + (3 J qdfi,fi G Mi*j = sup a/3 + H(a). 

The function H is concave, continuous on the interval defined earlier : 

q~ = lim inf (gW, 

q + = lim sup (g)r- 

In the case of a negatively curved manifold, this definition coincides with 

g~ = inf < / qdp,,p,E Aii 
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q + = sup I J q dp, p G Mi 

The function H is real analytic and strictly concave in ]g_,g + [. The function P is real 
analytic, strictly convex on R. Clearly, P(/3) > (3q+ for j3 > 0, and it is not very difficult to 

show that the limit lim^ >+oc P(f3)—f3q + exists and is nonnegativdl Similarly, P(/3) > (3q_ 

for j3 < 0, the limit linig >_oo — exists and is nonnegative. 

The pressure and the entropy appear naturally when studying large deviations for the 
Birkhoff averages of the function q. Denote Jt{p) the Jacobian of DG f going from E u (p) 
from E u (G t p). Define <p (the infinitesimal unstable Jacobian) by 

at |t=o 

On a manifold of dimension d and constant negative curvature —1, the function ip is 
constant, equal to d — 1. In general one can only say that it is a Holder function. The 
function ip is not necessarily positive, but it is cohomologous to a positive function, for 
instance (<p)t for T large enough. In what follows, we will assume without loss of generality 
that p> > 0. 

The two following large deviation results are due to Kifer [T3] . 



Theorem 3.1. [13], Prop 3.2. Let M be a compact manifold of negative sectional cur- 
vature. Let q be a smooth function on S*M . For T > 0, define the function (q)x = 
t I-t/2 $ ° G s d$ on S*M. Denote Li the Liouville measure on S*M. 
Then 

log f s , M e T ^dLi(p) 
lim — = Pr(q — (p). 

As a consequence, one also has : 

Theorem 3.2. [TB], Thm 3.4 (i) Let M be a compact manifold of negative sectional cur- 
vature. Let q be a smooth function on S*M. For T > 0, define the function (q)r = 
T I-t/2 1 ° G s ds on S*M. Denote Li the Liouville measure on S*M. 
Then, for any closed interval L C R, we have 

log Li {p e S*M, (q) T ( P ) e /} r f kA f \ 

hmsup — < sup < n KS {p) — / pdp,p GMi, / qdp G I > . 

T-^+oo T I J 2 J J 



For any open interval L C R, we have 
log Li {p G S*M, (q) T (p) G /} 

m int 

T 



log Li {p G S*M, (q) T (p) G /} r /• /" 1 

liminf 2 — > sup < h KS (p) — J ipdp,p G M.i, / qdp G / > . 

(mJ Let M be a compact manifold of dimension d, with constant sectional curvature — 1. 
TTien ^ can be rephrased as follows. Let q be a smooth function on S*M. For T > 0, 



This limit is equal to H(q + ). 
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define the function (q)x = \ J^t/2 1 ° G s ds on S*M. Denote Li the Liouville measure on 
S*M . Then, for any closed interval IcK, we have 

log Li {p E S*M, (q) T (p) E 1} 
lim sup — < sup {H(a), a E 1} — (d — 1) 

T >+oo T 



where H is the function H (a) = sup ^JiKsip): P E M.i,f q dp = aj . For any open interval 
I C Rj we have 

log Li {p E S*M, {q)r{p) E 1} 
hmmi — > sup {H{a), a E I } — (a — 1). 

The pressure and entropy functions also appear when counting closed geodesies 7 with 
a given g-average : 

Pr(g) = lim - log V e^ qd \ 

t — >+oo t — ' 

and as consequence 

lim sup - log (j 1 7, l 7 < t , I q dp^ E I \ < sup {H(a),a E 1} , 

t »+oo t I J ) 

(for a closed interval /), 

liminf - log (j {7, L < t, / q dp^ E I \ > sup {H(a), a E 1} , 
t I J ) 

(for an open interval I). See [T4] . 

In negative variable curvature, we will also need the following variant of Theorem 13.21 : 

Theorem 3.3. Let M be a compact manifold of negative sectional curvature. Let q be a 
smooth function on S*M. Let (ft be a smooth positive function. Denote Li the Liouville 



measure on 



S*M. For pE S*M andtE R, define T p {t) by ^ p{t) (j){G s p)ds = (d-l)t. For 



t > 0. define the function 



C-t/2) 



on S*M. 

Then, for any closed interval Id, we have 

log Li {p E S*M, (q)T(-t/2),T(t/2)(p) E /} 
lim sup 

t >+oo t 
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For any open interval I C E, we have 

log Lx {p e S*M, (q)r(~t/2),r(t/2)(p) e /} 

lim inf 

f >+oo t 

>( d _l )Bu p/^^_/^ >Ai6A<1) f gd/16 J 

{ J 4>dfi } (j)dp 2 J 

Proof. Define a flow (5 on S*M that has the same trajectories as G but different speed : 
G t (p) = G Tp<yt \p). For the new flow, the infinitesimal unstable Jacobian is equal to (d — 

1) j^^ - If A* i s an invariant probability measure of G, then dp = j^j^ is an invariant 

probability measure of G. Besides, their entropies are related by the Abramov formula : 

h K s{n) 



h K s(v) = (d-l) 



J (p dp 



where the entropies of p and p are computed with respect to G and G respectively. 

The theorem is then again an application of Theorem 3.4 in [13] for the Anosov flow 
G. ' □ 

4. Averaging 

We are now ready to start the proof of Theorem 11.41 We will work in dimension d and 
constant negative —1. The changes to make in order to get Theorem 11.71 are indicated in 
Remarks 14.31 and 16.11 

The following proposition is proved in |28j, §2 : 

Proposition 4.1. Let T > 0, there exists an invertible selfadjoint pseudodifferential oper- 
ator At E ^DO° such that 

A-\P + ihQ(z))A T = P + thOp h (q T (z)) + h 2 R T (z) 

for z E Q; with Rt E ^DO° depending holomorphically on z E VL, and with q T (z) E S 1 
depending holomorphically on z eVL, equal to (q)r — q + q z in a neighbourhood of p~ x Q) . 

The definition of our symbol classes S m and operator classes ^DO m is given in Section 

El 

We recall that the operator At constructed by Sjostrand is A = Op h (e 9T ), where 
g T = I [ ' - l) q o G s ds + ~ f (^ + l)qo G s ds 



2 Jo \T J 2 J_ T /2 V T 

on p~ l (|). The function qt solves {p ,9t} = q — (q)t- Exactly the same proof yields : 

Proposition 4.2. Assume M has constant curvature —1. Lete > andT = (1— 4e)| \ogh\. 
Define 5 = ^r. There exists an invertible selfadjoint pseudodifferential operator At E 
$>DO° 5 such that 

A T \P + ihQ(z))A T = P + thOp h (q T (z)) + h 2 R T (z) 
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for z G Q; with Rt G h~ 2S ^DO$ depending holomorphically on z, and with q T {z) G S} 
depending holomorphically on z, equal to (q)r — q + q z in a neighbourhood of p~ x (|). 

In what follows, we will restrict our attention to a region where \z — \\ = 0(h). As a 
consequence, we can write 

(4.1) P + ih Vh {q T {z)) + h 2 R T {z) = P + ih Ov h {q T ) + hR T (z) 

with q T = q T (|) = (q)r in a neighbourhood of p" 1 Q), and Rt(z) is a pseudodifferential 
operator depending holomorphically on z G Q, tending to zero when H — > and \z — \ \ = 
0(h). More precisely, 



R T (z) = (z - Q\z) + HR T (z) 



Rt G h 25 ^DO® depending holomorphically on z, and Q\z) G ^DO} depending holomor- 
phically on z. 

Remark 4.3. To treat the case of variable curvature, we should modify Proposition 14.21 
as follows. Fix a smooth function such that > (p. Define T p (-j), T p {— as in Theorem 
13.31 in a neighbourhood of p~ l (|). We have to choose <p smooth because we want T p to 
depend smoothly on p. In dimension d = 2, we have T p (^) G S$ with S = ^ (which may 
not be true for d > 2 since the unstable Jacobian no longer controls the derivatives of the 
geodesic flow). In Proposition 14.21 we now define A T = Op h (e 9T ) where 

on p" 1 (|). In the last sentence of Proposition 14.21 we replace (q)x by (q) T (-T) : T(^) + 
where r T = q-{p a , #T}-(g}r(-f ),r(|) satisfies r T G | log^| _1 5^ and r T } G | log ^l -1 ^. 
For d = 2, all the operators i? T etc stay in the same class as stated in Proposition 14.21 

5. Perturbations with controlled trace norm. 

In the following sections, we let z vary in a disc of radius 0{h) around |. We will write 
2z = 1 + £, ( = 0{h). We consider the operator (14. IK . that we write 

(5.1) VT = V T (z) = P + ihQ T + hRT(z), Q T = Op h {q T ). 

Note that we have 



in a neighbourhood of p G 1 (|). By Proposition 19.21 this implies 

(5.2) \\[P,Qt]u\\ <c(^ + 0(h 2 - 25 )^) \\u\\ + 0(h)\\(P - 

We now want to make a small perturbation VofV with a good control over the resolvent 
(V(z) — z)~ l , and over the trace class norm ||.||i ofV — V. 
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We construct a pseudodifferential operator Qt = Op h (q T ) G ^>DO] such that q T < q T on 
Po 1 (I) an d {Po,q T } = O [hj. In addition, we fix some e > and introduce an arbitrarily 
small 9 > 0, and we want q T = q T on p~ — e, | + e[) fl {q T < a — 39], and q T < a — 29 
everywhere on p~ — e, ^ + e[). For instance we can take q T = a(q T ) on p~ — e, ~ + e[) 
where a is real and smooth, a(E) < E, \a'\ < 1; a(.E') = -E if _E < a — 39, and a < a — 2$ 
everywhere. 

Remark 5.1. At this stage it is convenient to choose a positive quantization scheme Op^, 
in order to have Op ft (g T ) > Op^(g T ). 

Let < / G S{R), with / G C£°, where f(t) = J e itE f{E)dE is the Fourier transform. 
Put 

V = P + ihQ T + hR T {z), 

with 

(2P — 1 \ ~ (IP — 1' 

QT = Q T + f[ — (Qt - Q T )S ~ 



ft J^ 1 ^ ±,J V h 

The following proposition is proved in [28] for fixed T (and S = 0, that is, with standard 
symbol classes) . One can follow the proof of [28] line by line and check that it is still valid 
for T = (1 — 4e)| log ft |, e > very small : 

Proposition 5.2. Let P = — ft 2 ^-- Let Q = Q(z) G ^DO l have principal symbol q(z) 
depending holomorphically on z G VL, and be formally self-adjoint when z is real. Let 

V T = P + ihQ T + hR T (z) 1 Qt = Qt(1\ z = ^-^, ( = 0(h), 



be the operator defined in flOl . with Q T = Op h (q T ) G ^DO], and R T (z) G k^^DO^ + 
(z - \)^DO\. Let Q T = Op h (q T ) G ^DO], with q T = a(q T ) on p- l {\\ - e, \ + e[), where 
a is real and smooth, a(E) < E, \a!\ < 1; a(E) = E if E < a — 39, and a < a — 29. 
Put 

V T = P + ihQ T + hR T {z), 



with 



Then 



(2P — 1 \ ~ (2P — 1 

Qt = Qt + /( — —J {Qt-QtvI— — 



\Vt - Vt\\ < ft ( \\f\L sup {q T - f) + 0{h l ~ 25 ) 



and 

\\Vt-Vt\u <c d h 2 ~ d 



P(0) I (q T -q T )L,(dp) 

N-l f 
+ £ W/*(0)| / \D?(q T - f)\Li{dp) + 0(h N ^) 
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where Df k and D 2 p k are differential operators of degree < 2k, respectively on R and T*M. 

If we restrict z by assuming that for some continuous function a(E) > 0, defined on 
some bounded interval J containing 0, we have 



then for h small enough, (z — Vt) 1 exists, and we have 
Uz-Vt)) II < U 



h / n (®<o 



n ■ n 

The proof is identical to the proof in [28]. In Appendix [TOl we will give some details, for 
the reader's convenience. 

Corollary 2. Define 

Q h = ji - 2ch < fRe(z) < ~ + 2cfrj n {(a - 9)h < 5Sm(z) < 4^1100^} C C. 

For z G tin, the operator z — Vt is invertible, and 

\\{v T -z)- l \\< 



hi / 5t« 
en ' 

Corollary 3. For h~ small enough, we have 



< C^h^L, (V > a - 36} n {p; 1 Q) }) 



for small enough. 



We already defined 



6. Jensen's inequality 



n n = j~ - 2cft < Sie(z) < ~ + 2cfrj n {(a - 9)h < 5Sm(z) < 4:\\q\\ooh} C C. 
To finish the proof of Theorem 11.41 we also introduce the set 

Q h = I - - ch < Ue(z) < ^ + ch \ n {ah < Qm(z) < 3\\q\\och} C Q n . 
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For z G tin, we can write 

V T -z = (V t -z)(1 + K(z)) 

where K(z) is the trace class operator (Vt — z)~ 1 (Vt — Vt)- 

We can bound the number of eigenvalues of Vt in fin by the number of zeros of the 
holomorphic function g(z) = det(l + K(z)) in fl^. Let us call N(g,fln) this number of 
zeros. Introduce z = ~ + 2ih\\q\\ 0O . By the Jensen inequality [25] . 

(6-1) N(g,Q H ) <C (tog\\g\\ xAti -log\g(zo)\) , 

where the constant C does not depend on h (because the rectangles fin and fin can be 
transported, by translations and homotheties, to the fixed rectangles fli and fli). 
On the one hand, for all z e fin, 

\det(l + K(z))\ < exp\\K(z)\\i 

< exp (\\P T - z^WWVt - Pt\\i) 

< exp ( Cd J'1 ^d h \(d-l)~H(a-3e)-e](l-4e) 

~ \ Oh 

< exp (c M dh 1 ~ d h [{ ^ 1) ~ Hia ~™ ) ~ eK1 ~ 4e) ) ■ 

On the other hand, we also know that ||(1 + ( ) ) 1 1 1 < Cfc 1 : since z has 'large' 
imaginary part, Vt~Zq is invertible, and it is easy to get a bound \\(Vt — ) ~ 1 1 1 = 0{tT l ). 
We use the same calculation as in [28j and get 

\&e%{l + K(z Q ))- x \ = \&et{l-K{zv){l + K{z Q ))- l )\ 

< exp||K(^ )(l + ^o))" 1 ||i 

< expii^iMKi + ^o))- 1 !! 

< exp (c , /i9ie ,^ 1 -^[( d - 1 )- H (^ 3e ^ £ ]( 1 - 4e )) 

so that 

| det(l + K{z )\ > exp (^-C ftqfi4 h l - d h^ d ' 1) - H{a ^ e) - t]{l - Ae) ^ . 
This, combined to (|6.1I) . yields 

\N(g,fl h )\ < Ch l-d H l(d-l)-H(a-28)-e]il-ie) 

Since 6 and e > are arbitrary, we have proved Theorem II .41 

Remark 6.1. Starting from Remark 14.31 the proof of Theorem 11.71 goes exactly along the 
same lines. We find 











|log 


h\ 



limsup M — — t — < H(a) 

h — >o 



where H(a) = (d—l) sup i ^yfj^ — + 1, /i € M. i , / qdfi = a\ and is as in Remark 
14.31 Letting <fr converge to <p uniformly, we obtain Theorem 11.71 
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7. About Question 3 

In this section, we consider a particular case of the problem d 1 . 5 f) in which the trace 
formula is exact. We then try to investigate Question 3 on this example. 

Let M be a compact hyperbolic surface : M can be written as M = T\H, where H is 
the hyperbolic disc and T is a discrete subgroup of the group of hyperbolic isometries. Let 
[uj] £ if 1 (M, C) be represented by the harmonic complex valued 1-form uj. Introduce the 
twisted laplacian 

A w / = A/-2(w,#)+|M£/. 

Studying the large eigenvalues of A w amounts to studying a fixed spectral window for the 
semiclassical twisted laplacian 

-H 2 ^ = -h 2 ^ + h 2 (u, d.) - h 2 ^, h^o. 

The "usual" selfadjoint case is when uj has coefficient in iR. We shall instead be interested 
in the case when uj has coefficients in R. The operator falls exactly into the case studied 
in §1.21 with q(x,£) = ((jj x ,£). The geodesic flow is ergodic, and Sjostrand's result tells us 
that "most" eigenvalues of —h 2 ^- such that ^te(z) £ [| — CH, ^ + Ch] have imaginary part 
3m(z) = o[h). Equivalently, "most" eigenvalues of -A w such that 9fte(z) £ [A — CyX, X + 
Cy/X] have imaginary part 5sm(z) = o(vA). 

We rephrase Question 3 as 
(Q3') Uu^O, is it possible to have — ► as K — ► and $le(z) £ [§ - CH, \ + CH], 

z £ s P (-h 2 ^ ? 



h 

.\ ? 

2 



Conjecture : I conjecture the opposite : if uj ^ 0, then there is a sequence h n — ► 0, 

[I - Ch n , I + CH n ] and 



z n £ Sp(-h 2 n ^) with Me(z n ) £ [f - Ch n , \ + Ch n ] and 0. 



As is usual in hyperbolic spectral theory, we introduce the spectral parameter r : if Xj is 
an eigenvalue of — A U) we denote Xj = \ +r 2 . Yet another way to phrase Question 3 is to 
ask whether it is possible to have Qm(rj) — > as §fte(rj) — >■ oo. Sjostrand's results say 
that Sm(r 3 -) is bounded and that Qm(rj) — > for a subsequence of density one. But I 
naturally believe that it is impossible to have Qm(rj) — > for the whole sequence, unless 
uj = 0. 

Recall the Selberg trace formula [27], valid for uj £ f> rl (M, iR) : 



(7.1) = P/>)rtanh(Hrfr + ^^f /(/ 7 ) 



7/' 

2 



for any function / on IR, even, smooth enough, and decaying faster than any exponential. 
On the right hand side, the sum runs over the set of closed geodesies (equivalently, the 
set of conjugacy classes in T). If 7 is a periodic geodesic, we denote / 7 > its length, or 
period; and l Jo is its shortest period. 
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Proposition 7.1. The trace formula holds, under the same assumptions on f , if u G 
H 1 (M, R) . 

Recall that the Fourier transform is defined by /(r) = J e vru f(u)du. 

Proof. Take u G i/ 1 (M, R). We consider the operator A ZU} for z 6 C. The argument of 
Section [2] shows that this operator has discrete spectrum (eigenvalues). 
The right hand side of the trace formula reads 

, , Area(M) f + °° P . . . , ™ e *A w L 

(7.2) 7^/ /(r)rtanh(vrr)rfr + V— -^/(Z 7 ) 

47T <—? smh - 



31 

7 2 



and clearly is an entire function of z. 

The left hand side is f( r j{ z ))- To check that it is an entire function of z, we first 
note that f(r), being an even entire function, can be written as g{\ + r 2 ) where g is entire. 

Thus f(rj(z)) = g(Xj(z)), where the Xj(z) are the eigenvalues of — A ZUJ . 

If z is restricted to a bounded subset Q of C, we note that the Qm(rj(z)) are uniformly 
bounded. To that end, we write — Xj — \ + r| with r n = x + iy. The eigenvalue equation 

- A zw / = A 3 -/ 

with / normalized in L 2 implies both equations 

i + * 2 -y 2 = | lV/f + 2/^ J {uJ ,df)f + (f3 2 -a 2 ) J\\u\\l\f\ 2 

and 

2xi/ = -2m J(u,df)f-2af3 j ' \\u\\ 2 x \f\ 2 

if we decompose z = a + GR + iR and = x + G R + zR. If a and (3 stay bounded, 
it also follows that y must stay bounded. 
Besides 

A~ 2 |j {n, < Ue{r n {z)) < X} 
is bounded uniformly for A > 1 and z staying in a compact set of C (the arguments of 
[28] . §4, or of our Sections 01 El El are locally uniform in z). Since / is rapidly decreasing 
in each horizontal strip, it follows that the sum J2j f( r j( z )) i s the uniform limit of the 
partial sums Yli\Ue{rj(z))\<\fi r A z ))- But f° r a gi ven \ this is a holomorphic function of z 
(in the open set {z, $te(rj(z)) ^ X for all j}), since A zu) depends holomorphically on z and 
Yli\®e{rj(z))\<\ f( r j( z )) can be defined by holomorphic functional calculus. 

This shows that Ylj f( r j{ z )) ^ s a ^ so an entire function. Both sides of (17.21) coincide for 
z G iR (by the usual trace formula), thus they must coincide everywhere. □ 

Introduce some parameters a,R,T> 0, and take 



/(«) 



27TCT 



(m-T) 2 . (u+T) 2 . 

e — 5^e mi? + e — 5^e~ mi? 
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SO that 



f(r) = e-*-^ e~ lTr + e~^ r+R ^ e lTr . 

This yields : 

(7.3) e-^- R)2 e- lTr > + e -£to+*) s 



Area(M) 
Air 



/+OO 2 2 

rtanhvrr I'^'^e^ + e -¥('+*>V Tr <ir 
-oo 



eA"^ 1 



sinh ^ \phxo 

7 2 



" 2^ e 7 + e e 7 



We want to bound from below the right hand side. We hope that this bound will tell us 
that exp(±iTr\,) cannot be too small on the left hand side, giving some information on the 
imaginary part of rj, for 9fte(rj) ~ R. On the right, the idea is that the main contribution 
should come from the geodesies with Z 7 ~ T. We want to choose R so as not to be bothered 
by the oscillatory terms e ±tllR . For that purpose, R and T will be related in the following 
manner : 

Lemma 7.2. [21], Lemma 3.3, [HJ. 

For any M > 1, there exists R G [M, Mexp(exp(5T))] such that cos(Rl 7 ) > \ for every 
closed geodesic 7 with / 7 < 5T. 

In the sequel we take M = exp(exp(cT)), (c > arbitrary) to ensure that T is of order 
log log R. We note that this relation between R and T is independent of u>. This will allow 
us to modify slightly our initial problem by extending it to the case where u can depend 
on R (or T). More precisely, we want to consider the case when uj = Q(R)uj , where uj q is 
fixed and Q(R) > 1 is allowed to go to infinity with R at a reasonable rate. 

This means that we consider a slight generalization of (I1.5P (the motivation should 
become clearer in §T.2j) : 



7.1. A more general problem. We consider a spectral problem of the form 

(7.4) (V -z)u = 

where 

V = V(z) = P + iH0(h)Q(z), P = -^- 

where z6fl = e i] - flo ' So[ ]E min , E max [, with < E min < \ < £ max < +00, < s < f . We 

will assume that Q(z) G ^/DO 1 depends holomorphically on z G Q, and that 8(h) is some 

real valued function such that 8(h) > 1 and K6(K) — >■ 0. We have in mind 6(h) = \ log(^)|. 

h — >o 

Finally, we assume that Q is formally self-adjoint for z real. Again, we call S the "spectrum" 
the set of z for which the equation (V(z) — z)u = has a solution. 
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The results of §1.21 can be generalized as follows : for any E min < Ei < E 2 < E max , 

1 



(7.5) {t {z G E, E 1 < 3te(z) < E 2 } 



(2nh) 



d 



dxd£ + <D(h9(h)) 



One can show that has to stay bounded for z G E. Taking #(/?,) = | log(ft)|, £7 

-E — ch9(h) and E 2 = E + cH9(h), one has 



for z G E such that £i < 3fte(z) < E 2 . Assuming that the geodesic flow is ergodic on 

p~ 1 {E}, and for 9(h) = \ log(ti)\, then for any e > 0, any c > 0, 

(7.6) 



z G E, E - cft0(ft) < Vte(z) <E + cH9(H), £ [q E -e,q E + e]j = 9(h)o(h 1 - d ). 

Remark 7.3. The paper [28] only treats the case 9(h) = 1. But the method of [28], §5 can 
be adapted in a straighforward way to show the following : consider the spectral problem 

(7.7) (V -z)u = 

where 

h 2 A 

V = V(z) = P + ih9Q(z), P = — 

where z G Q = e^- s °' s ^}E min , E max [, with < E min < \ < E max < +oo, < s < ~ . Fix 
some e > 0. Then 

1 



(7.8) ft {z G E, Ex < Me(z) < E 2 } 



(2-nhY 



dxd£ + 0(h9) 



uniformly in all 9 > 1 such that 9h < e and E\, E 2 such that E min < Ei—2e, E 2 +2e < -E max , 
\E 2 -Ei\ > M. . 

For (17.61) (and 9(h) = \ \ogh\), the generalization of the proof in [28] is without surprise, 
but requires some rather technical changes : we will not prove it here, but still feel allowed 
to ask about Question 3 in this generalized setting. We note that to extend (|7.6I) to more 
general 9(h), some analyticity assumptions would be required, as in [9]. 

We will focus our attention on the operator 

_ h2 = _ h2 A + d) _ ft 2 9(ft)2 W ) h _^ Qj 

when u has coefficients in R. 

We generalize Question 3 as 
(Q3") If uj 0, prove that there exists a sequence h n — > 0, and z n G Sp(—h^ 0( ^ n) " ) with 
Ke(z„) G [i - Ch n 9(hn), \ + CMW], such that g[g} 0. 
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7.2. Heuristic discussion of the parameters a, R, T. We start again from the trace 
formula (17.31) . considering the case where uj = Q(R)u , Q(R) = 9{R~ l ) > 1. Here R^ 1 is 
going to play the role of the small parameter H. The form uj is fixed, and we normalize it 
to have stable norm \\uj \\ s = 1. 



(7.9) W^Ws = sup < / uj dfi, /i E M. i 

(7.10) = sup^— . 

The first line can be considered as a definition of the stable norm (valid for a general 
compact riemannian manifold M), whereas the second line holds on negatively curved 
manifolds because of the density of the closed geodesies. Using the definition Vy{uj) = 

sup foiKsin) + Js*m uj ' dpi ^ e -^f }' ^ * s n °t difficult to show that 

lim Pr(ta>) — |t|||u;|| s = sup < h KS (n), \i E Aii, / udji — ||w|| a f . 

I 2 J5*Af J 

On a surface, the right-hand side vanishes Besides, for any T one can find a closed 
geodesic 7 with L t E [T — 1, T], and such that 

(7.11) 4^ > n^mi + 0r( i)) = (1 + 0r (i) ); 

where 0^(1) goes to as T approaches +00. Simply recall that for any < 5 < 1, 
] ^logtl{7,/ 7 G[T-l,T],^>(l-5)} 

where 



(7.12) lim 25 — = (1 - 5) > 0, 



if(a) = sup < hKs(fi) ,y>E Mi, I u dfi 

I 2 JS*M 



a 



The function if is continuous, concave on [—1, 1], real-analytic on ] — 1, 1[ [4]. And again, 
H(—l) = H(l) = on a compact surface [2]. 

In (I7.3p . we have not said yet how a will depend on R and T. For the moment, let us 
decide a priori that a should be such that the term 



Area(M) 
An 



r°° •- 2 2 

^(r-R) 2 e -iTr + ^^{r+R) 2 ^Tr 



r tanh nr 



dr 



is negligible compared to the sum ^ . Remember that i? and T are chosen so as to satisfy 
Lemma 1731 Then, fixing 1 > 5 > 0, the right hand side of (17.31) should be bounded from 
below by 

(7.13) a -l e TH(l-6)-T/2 e e(R)(l-6)T e -^ 
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We want to use the fact that this grows quite fast with T. On the other hand, looking at 
the left hand side of (17.31) . we cannot hope to do better than to bound it from above by 



(7.14) 



{j, \^e{rj) -R\< a' 1 } e V su P, Sm (^)V sup ^ m(r ^ 



where each time the sup^ should be restricted to the indices j such that |3te(r 3 -) — R\ < a 

This heuristic argument would give an inequality 
(7.15) 



-i 



{j,\3te( rj )-R\ <a- 1 }e £ 



■ sup Sm(rj) 2 T sup • |5Jm(rj 



> a 



-l p TH(l-5)-T/2 p Q{R)(l-5)T 



e 2 



obtained by comparing the lower bound (17.131) and the upper bound (IT. 14j) . Again, the 
hope is to compare the powers of e T on both sides to prove that sup^ |9m(Vj)| cannot be 
arbitrarily small. 

Consider the case Q(R) = 1, which is the case we were originally interested in. If we 
take a to be a constant, then by Weyl's law we have (j {j, |9fte(rj) — R\ < o" -1 } ~ R > 
exp(exp(cT)). In this case (17.151) cannot bring any useful information. On the other 
hand, if we want to choose a such that |j {j, |3?e(r 3 ) — R\ < cr^ 1 } is bounded, we are led to 

^ 2 

take a ~ R; in this case the term sup J 9fm ( r ^*) will be too large to yield any interesting 
information. 

We see that the method only has a chance to work if TQ(R) ^> logi?. From now on we 
take Q(R) > logR, and always such that R" 1 Q(R) — > 0. We also take a~ 2 = CQ{R) with 
C large. We must note that the parameters rj correspond to the eigenvalues of — Ae^, 
and thus they also depend on R. 

7.3. Proof of Theorem 11.121 The right hand side of (17.31) is easy to understand. The 
term 

" + 0O 



(7.16) 



r tanh nr 



(r+R) 2 e iTr 



dr 



is 0(a l R), whereas the ^ has modulus greater than 



(7.17) 
1 
2 



f 



<5T Sinh i 



2ttct 



2a 1 -\- e 2CT 2 



f U) J 

+ — ■■ 



sinh 



llXO 



(h-T) (h+T) 
2a 1 -\- g 2a 2 



> 



1 



; j t I, 



lo 



2 T-t^l<T Sinh 2 ^ 27r(7 



ih- T r 



+ e 



eh l 7o 



> 



7 >st sinh i V 2vra 



e 2<r 2 



+ e" 



cos(/ 7 i?) 

cos(Z 7 i?) 



e TH(l-S)-T/2 e e(I{)(l-5)T e -^ a 



0(1) 



(using (17.121) ). and thus is much greater than the integral (17.161) . To get the last 0(1) we 
have used the exponential growth of the number of closed geodesies. The left hand side of 
(17.31) is more complicated to bound from above, since the rj now depend on R. 
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Proposition 7.4. Take Q(R) > log R, and such that R~ 1 Q(R) — ► 0. Take a' 2 = CQ{R) . 
Let f(R) be such that a 2 f{Rf > TQ(R). Then 

|^ e -4fo-*)V iTr i| < l&efo) -R\< /(#)} e 4 s ^^i)V su P^ m fo)l + 0(1), 

where the sup- are taken over the set of indices j such that \dte(rj) — R\ < f(R). 
Proof. 

Lemma 7.5. VFe /iai>e an a priori bound \^sm(rj)\ < cQ(R), where c depends only on uo Q . 
Indeed, let rj = x + iy and / G L 2 (M) be such that \\f\\L 2 — 1 and 

A 

Taking the scalar product with /, we get 

(7-18) ^ + * 2 -y 2 = /|V/| 2 - Q{R? J \\uj \\ 2 x \f\ 2 

and 

(7.19) 2xy = -2Q(R)i J(u ,df)f. 



A / + 2e( J R)( Wo , df) - e(R) 2 \\u \\ 2 x f =(- + r 2 )f. 



Equation (17151) yields < cQ(R)J f |V/| 2 . Equation (ITlSll implies that x 2 > f |V/p 



c 2 6(i?) 2 . If |x| > |y / | V/| 2 then we are done, by (17191 . If |x| < |y / |V/| 2 , then fl7l8l) 
implies that / |V/| 2 < 2c 2 ®(R) 2 and that y 2 < \ + 5c 2 e(i?) 2 . The lemma follows. 

We now break the sum Ylj e ' r ^ rj ~ R ^ e~ tTrj into three parts : / = Y^j^ste{r j )<R-f{R)^ 

11 = 22j,\^e( rj )-R\<f(R) and 111 = J2j,Ue{ rj )>R+f(R)- 

The last sum III is bounded by 

e 4c 2 e(R) 2 e cre(R) e -^{^i)-R)\ 

j,3te( rj )>R+f(R) 

We decompose this sum into J2 n >o T,R+f(R)+n<^e(r j )<R+f(R)+n+x> and b y We y rs law in the 
form ( 17.81) . this is dominated by 

2 P+OO 

e cTe(R) J-(R + n + f(R))Q(R)e-^ n+f ^ 2 < e cT6 ^Q(R) / (R + 



x)e 2 dx 

my 



and with our relations between T, R, f(R) and a, this last quantity is 0(1). 
Concerning the first sum /, we bound it by 



(7.20) e ^e(R)* e cTe(R) e 2 

jMe(rj)<R-f(R) 
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The subsum ^ m e (r)<-R+f(R) can be treated as above and shown to be 0(1) (using Weyl's 
law in the form (17.81) ). and we only need to concentrate on Ylj,\9ie(r-)\<R-f(R)- To bound 
this sum, we first need a control the number of terms. 

Lemma 7.6. {j, \^e( rj )\ < R} = 0(R 2 Q(R)). 

To that end, we use again the trace formula and write : 

_J Area(M) f + °° , _ , 2 , ^e e(fi) ^°L R _ r2/2/2 

\ e 2h2" = 1 — '- I rtanh(7rr).e ^dr + > ; — -^^e t ; . 

^ 47r ,/_ 0O ^— ' sinh^- V27T 

J T * 

On the right, the term X] 7 is clearly o(l) whereas the J term is of order i? 2 . On the left, 
we break the sum into ^ |gte( r -)|<R anc ^ \9te(r )\>R- ^ s above, we can use Weyl's law in 

, , r | (fl+n) 2 

the form (jZSJ to show that the Ei,|SR e (r 3 )|>ii e ^ is °J2 n >o( R + n + 1)©(-R)e = 
0(i? 2 9(ii!)). Thus we have 

e 2^ 1 e 1 W— L - = 0{R 2 Q{R)). 

j,We( rj )\<R 

Remember that l^mfo)! < cQ(R) and that R^ 1 Q(R) — ► 0. Thus, for |3te(rj)| < R we 
can write e ^ = 1 + O(i2 -1 0(i2)). This yields 

, »e(r,) 2 -am(r,) 2 

e-^O'J^r,)! < ^Xl + O^e^))) < e (l + OiR^eiR))) 

j,\3terj\<R 

= 0{R 2 Q{R)) 

and finishes the proof of Lemma I7T61 

Now, we go back to the sum Yljme(r )\<R-f{R) m ^7.201) . and we see that it is bounded 

by 

R 2 Q(R)e^ c2@ ^ 2 e cT0 ^e^ 2f ^ 2 / 2 = 0(1). 
This ends the proof of Proposition 17.41 □ 

We can now come back to (17.31) . Noting that a 2 supj ;$$m(rj) 2 = 0(1), Proposition 17.41 
shows that the left-hand side of the trace formula ( 17.31) is bounded from above by 

Cfl {j, \®e{ rj ) -R\< f(R)} e Tsup ^ |3m(r ^ + 0(1) < CRf(R)e(R)e Tsup i |3m{r ^ 

where the sup is taken over all j such that |9£e(rj) — R\ < f(R), and where we have used 
again (EH}. We can of course assume, without loss of generality, that f(R) = 0(R). 

On the right hand side of (17.31) . the J is 0(i?9(i?) 1//2 ), and the ^ 7 is bounded from 
below by 

1 _ e TH(l-S)-T/2 e e(R)(l-S)T e -^j 



87T(T 
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as in ( 17.171) . Writing 

CRf(R)Q(R)e Ts ^^ m ^ > a -^ e TH(is)-T/2 e e(R)(i-5)T e -^ ? 

remembering that Q(R) > \og(R), Q(R) = o(R), a~ 2 = cQ(R) and T x log log R, we see 
that necessarily 

suplSmfo)! > (1-25)B(R). 



3 



This finishes the proof of Theorem 11.121 



8. The arithmetic case. 

Let p > 3 be a prime, p = 1 (mod 4), and A > 1 be a quadratic non-residue modulo p. 
We set 

r = T(A, P ) = {( v ° + mV %- y2 ^ + v '^),yo,m,y„y,£z). 

I V 2/2 VP - 2/3 V^P Z/o - ViVA J J 



It is a discrete cocompact subgroup of SL(2, R) which contains only hyperbolic transfor- 
mations [8]. We consider the hyperbolic surface M = T\M. In M, the lengths of the closed 
geodesies are the logx m , where 



x m = 2m 2 — 1 + 2m V m? — 1, m G N. 

We define 

7, l~f=\ogx m 

We now follow very closely the approach of [8], pp. 304-314. We introduce an even function 
k on R, whose Fourier transform is nonnegative and compactly supported in [—1,1]; we 
also assume that k > 1 on [— ~, -]. We define 

K Q (r) = A;(r)[e w + e~ iar }. 

We write again the trace formula : for all t > 0, 

(8.1) ^ifafo - t) + K a { rj +t) = Are ^ M ^ I rtanh(7rr )[K a (r -t) + K a (r + t)]dr 

+ 2 V ^h*K a {l 7 ) cos(« 7 ). 
- sinh -J 

We will bound from below the right-hand side (averaged in t) to obtain information on the 
left-hand side. Denote 



Sa{t) = Y ~T~~T^o l {lrf) cos(t/ 7 ) = 2 ^ l/ ^ m \ 1/2 k a (\ogx m ) cos(tlogx r 

7 Sinll 2 e a - 1 <x m <e Q + 1 Xm + Xri 



DAMPED WAVE EQUATION 27 

Proposition 8.1. Let a = 2(3 log T — C, < (3 < 1 and C large enough. Then, 

Although we are really interested in the quantity J 2T _ Tf3 \S a (t)\ 2 dt, the reason for intro- 
ducing the regularizing factor ^1 — j is exactly the same as in [8], p. 315. 

Proof. Introduce the notations 

rv> ' I ry-> » 

z/( m ) = /i(m)i^ a (logx m ). 
Divide the integral J = J^T^ ^1 — ^-ff^ J |£q(0I 2 ^ into / = 7i + -Z2, where 

h= J~] v( m f / 1 ^ — cos 2 (tlogx m )dt 

and 

Ii = 2 r](m)r](k) (l ^ — J cos(tlogx m ) cos(t logx fc )dt. 

The idea is that the x m ,Xk, with x m 7^ x^, are well-spaced, implying that the oscillatory 
integral J 2 is small compared to I\. 

Lemma 8.2. For T > 1 and AeR, 

-2T+T' 3 



12T-T1 3 \ 1 
Let us first consider ii. 

*2T+TP / \ + oti\ r 2T + T>3 



mm 



y/3 



1 



f^A |t-2r| ^ 2 / 2T+TP / ^-2T| \ l + cos(2tlogx r . 



T /3 

T + 0(T^loga ; - 2 ). 



Hence, 



(8.2) 7i > ^ + (1)^ Yl 

(T p , A \- z/(m 



(8.3) 

(8.4) > Cl T^- a 



1 ,„,<e Q + 1 + Xm 

2 



2 

V2N2 



e a - 1 <x m <e a + 1 
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The right-hand side of (18.41 ) will be estimated later. We now turn to I2, and want to show 
that it is much smaller than I\. The integral 

J ^ (l -^p — J cos(tlogx m ) cos{t\ogx k )dt 

is smaller than t^tt, — ; ™. We can ensure that 



\logx m - \ogx k \ 

for e a_1 < x k < x m < e a+1 , by choosing a in an appropriate range : writing 

•Em 2 4 



logx m - logx m _i 



-) 



we see we have to take a < 2(3 log T — C (C large). More generally, we have by the 
intermediate value theorem 

I \ogx m — logXfc| > Ce~ a ' 2 \m — k\. 
The analysis done by [8], pp. 310-311, can be applied verbatim to show that 

(8-5) \h\<^ E ^ 

e a - 1 <x m <e a+1 

Comparing (18.41) and (18.51) . we see that 

< ± 

1 21 - 100 

provided a < 2/3 log T — C with C sufficiently large. Thus, 

r 99 r 
I > h. 

- 100 

To complete our estimate for /, we must return to equation ( 18.41) . Clearly, 

E u ( m ) 2 ^ E fi{m) 2 . 

e a - 1 <x m <e a + 1 e a-l/2< Xm < e a+l/2 

We write the Cauchy-Schwarz inequality, 



E ^( m ) 

ea -l/2 <x <£a +l/2 



2 



< I V ! V Mm) 2 

' 1 / 2 <x m <e a + 1 / 2 I \e a - 1 / 2 <x m <e<*+ 1 / 2 



But 

E M(™)= E eW*>Ce« p *'>, 

ea -i/2< a . m < ea +i/2 7 ,a-l/2<i 7 <a+l/2 

see [20], p. 117. On the other hand, 

E i = ^ /2 )- 



- 1 / 2 <x rn <e a + 1 / 2 
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We obtain this way 

e a-l/2< x . m < ea +l/2 

We have proved 



2T+T /3 , . 



□ 



This implies 

for some t e{2T -T?,2T + 

Now consider the integral J r tanh(7rr)lf a (r — £)dr or J r tanh(7rr)ii" Q ,(r + i)dr in (18.11) . 
We write 

rtanh(vrr)ir Q (r - t)dr = [ r tanh(vrr)A;(r - t)[e ia(r - 4) + e^^dr. 



To evaluate f(r + t) tanh(7r(r + t))k{r)e loa dr ', we shift the integral over R to an integral 
over (| — e)i + R, and we find that the integral is (9(te~' a '^~ e )) for any e > 0. 

Lemma 8.3. J r ta,nh(irr)K a (r — t)dr or J r t&nh.(nr) K a (r + t)dr = 0(te~\ a ^~ e ^) = 
C(t 1+e_/3 ) for anye>0. 

We finally turn to Ka( r j ~ t) + K a (rj + t). Fixing a small e > 0, one sees using 
Weyl's law, the fact that k is rapidly decreasing in any horizontal strip - and the fact that 
the $sm(rj) are bounded - that 

j j,\9le(rj)-t\<t e 

Similarly, 

Y,K a {r 3 +t)= J2 K a ( rj + t) + 0(t~°°). 

j 3,\toe(rj)+t\<t* 

We see that 

\Y^K a {rj -t)+ K a (rj +t)\< Ce asup J l 9 ™Mt 1+e < Ct 2f}aup i l 9 ™^V+ e 
j 

where the sup ■ is taken over the j such that |9fte(rj) ±t\ < t e . 

We have proved that there exists t G [T — T^,T + T' 3 ], and tj with |9fte(r,-) — 1\ < t e , 
such that 

In particular, if T is large enough, this implies 

sup {|9mfa)|, IStefa) ± T| < T^} > Pr(w) - ^ - 
and this proves Theorem 11.91 
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9. Symbol classes 
Following [6], for any 0<<5<l/2we introduce the symbol class 

(9.1) = {a G C°°{T*M), \d^a\ < C a3 h~ & \ a+ ^ (£)"H/?I j . 

If 5 = we will just denote S m . We will denote Op h (a(x, £)) = Op(a(x,h£,)), where Op 
is a quantization procedure on M. The quantization of any a G S$ leads to a bounded 
operator on L 2 (M) (the norm being bounded uniformly in h), see [6j. We will denote 
*DO? = Op h (S?). 
We use : 

Proposition 9.1. Let a G Sf,be S% with < 6' < 5 < 1/2. Then 

(i) Op h (a) Vh (b) - Vh (ab) G h 1 - 5 - 5 ' Op ft (ST +n - 1 ). 

(ii) [Op ft (a), Op R (6)] - f Pa ({a, b}) G h^-*) Vh {Sf^- 2 ). 

We also use a local form of the Calderon-Vaillancourt estimate [6] : 

Proposition 9.2. There exists K G N depending only on the dimension of M , such that 
the following holds. Take A = Op h (a) where a G ^DOj. Let I be an open interval ofM. + 
and let A belong to I. Then, there exists C > 0, and C(a, A) depending on a finite number 
of seminorms of a (uniform in A if it stays inside a compact subset of I), such that, for all 
u ■ L-[.\I). 



\Au\\ L 2 < C ( sup \a 



+ J2 hk SU P \D 2k a\] N| L 2 + (7(a,A)||(P-A)u|| L 2. 

k=l Po\l) J 



In fact C(a, A) is controlled by the supremum norm of — ^ and a finite number of its 
derivatives outside p^ 1 (I) . Similarly we have 

\{u,Au)\ <C ( sup \a\ + J2 hk SU P \D 2k a\\ \\u\\ L 2 + C(a, X)\\(P - X)u\\ 2 L2 . 

10. Appendix : Sjostrand's proof of Proposition 15.21 

In order to prove Proposition 15.21 we first want to bound the norm and trace norm of 

/2P-1\ ~ /2P-1 
/(— E— ) Ot-Ot)/' 



h J^ 1 ^ ±,J V h 

We write a Calderon-Vaillancourt type estimate, 



\\(Qt ~ Qt)u\\ < I sup (f - f) + 0(h^ 25 ) ) \\u\\ + 0(1)\\(2P - l)ul 



DAMPED WAVE EQUATION 31 

where 5 is as in Proposition O Besides, \\(2P - l)f ( 2 ^ zl ) \\ = 0(h). It follows that 



/ 



'2P 



h 



{Qt — Qr)f 



'2P 



h 



< 



sup (q 



f) 



+ 0(h 



l-28\ 



For the trace class norm, we need to be even more careful than in [28]. Instead of using 
the Garding inequality, we use the existence of a positive quantization. If we choose such, 
we have directly that / ( 2P J~ X ) (Qt — Qt)} — m the operator sense. Thus, 



/ I ^"J— ^ (Qt - Qr)f 1 ~ P ~ ' 



h 



h 



^{^(Qt-Qt^ 

Trf(^^\ (Qt-Qt) 

Tr^- / P{t)e lt ^{Q T -Q T )dt. 
In J 



Writing the expansion of e ~ as a Fourier Integral Operator, writing the trace as the in- 
tegral of the kernel, and applying the stationary phase method in the time-energy variables, 
we obtain an asymptotic expansion 



Tr 



2n 



C d h 2 ~ d 



P(t)e u ^(Q T -Q T )dt 



r JV " 1 r 

P(0) / (q T - f)LUdp) + £ K k D?P(0) / Df(q T - f)L,(dp) + 0(h N ^ 



-25) n 



where D 2 and D 2 are differential operators of degree < 2k, respectively on R and T*M. 

Note that the term h k D 2k } 2 {0) J^/n Df{q T - q T )Li(dp) is a 0{h k ^- 25 ^)Li (q T ^ q T ) = 

o(l)Li(g T 7^ q T ). This proves, in particular, CorollaryEl 

To finish the proof of Proposition 15.21 there remains to study the invertibility of z — Vt- 
Recall the identity 



\\(A + iB)u\\ 2 = \\Au\\ 2 + \\Bu\\ 2 + i(u, [A,B]u), 
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if A, B are bounded self-adjoint operators. Thus, 
2\\(P T - z)u\\ 2 > \\{P + thQ T - z)u\\ 2 - 0{h 4 ^){\\{2P - l)u\\ 2 + \\u\\ 2 ) 
> || (P - fte(z))u\\ 2 + h 2 \\ (^HM _ Q T ^j + [ Pj q t]m ) 
-C(^ 1 - <5 ))(||(2P-l)n|| 2 + ||?i|| 2 ) 
= || (P - Ke{z))u\\ 2 + h 2 \\ f^p- - Q T ) n|| 2 

+ (o(l)^(l + ||/|| 2 J + 0(^ 3 - 25 )) INI 2 + 0(h 2 )\\(P - Ke(z))u\\ 2 
We have used (15.21) (or Proposition I9.2p . and the same for Qt- We find that 

(10.1) VS\\(V T - z)u\\ > ||(P - ®e(z))u\\ - (0(A) + 0(ftH)||u||. 
On the other hand, we have 

(10.2) $Sml-{z-V T )u,u 



^sm(z) 



h 



Q T )u,u) + Oih 1 - 25 )^ + \\($e(z) - PH|)||w| 



$Sm(z) n t ( 2P-\ \ ( ~ w/ gP-l . t 

Qt + / — r — (Qt - Qt)/ — = — ) ) u, u 



h \ h J \ h 

+ 0(h 1 - 2S ){\\u\\ + \\(Me(z)-P)u\\)\\u\ 

-r Qt + /( r I (Qt-Qt)|w,w 

+ 0(n 1 - 25 )(|| M || + ||(Ke(z)-P)n||)|| M || 

P - »e(z) 



2 sup (g T -g T )||/||oo||/ / ||oo + ^ 1 - 2!i 



it 



/l 



-It 



by the same trick as in [28], (3.19). Recall that we are interested in a region where 
z-\ = <D(K). 

Let a{E) > be a continuous function defined on a bounded interval J containing 0, 
and restrict z by assuming that 
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near p Q 1 (|) , £ J (where ( = 2z — 1). It follows from the Garding inequality that for 
such z 



Qr + f 



/23?e(; 



> I a 



(Qt ~ Qt) u, u 



0{h 



l-2(5^ 



\u\\ 2 - 0{l)\\u\\\\{P -$ie{z))u\ 



Using this in (1 1 . 2 [) . we get 



(10.3) 3m ( ^(z-V T )u,u ) > ( a 



- 2 sup (9 i -r)ll/||oo||/ , ||oo + 0(ft 

Reasoning as in [28], (3.25) and (3.26), we find finally 
1 



l-2<5\ 



U 



P - 



h 



-u 



h 



z — Vt)u, uj > ya 
->/3(2 sup (q T - f)\\fU\f\ 



1-25^ 



■ft 



l-25\ 



U 



h 



-II 



(0( 



T 



u 



and 



T 



0{h 



it 



< 



1 + 2^ sup (<? T -9 T )||/|Uin 



l-2(5> 



V T -z 



h 



-u 



which finishes the proof of Proposition 15.21 
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